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ABSTRACT 


The  problem  of  determination  of  the  transverse 
deflections  of  plates  under  conditions  of  pure  bending  in 
the  longitudinal  direction  is  considered  for  plates  with 
a  bi-trapezoidal  transverse  cross  section.  Previous  in¬ 
vestigators  have  considered  the  transverse  deflections  of 
plates  with  rectangular  and  double-wedge  shaped  cross 
sections.  This  thesis  extends  this  work  to  include  the 
case  of  a  bi-trapezoidal  cross  section.  With  a  suitable 
change  of  parameters,  the  developed  theory  is  shown  to 
approach  both  these  known  results. 

Experiments  carried  out  on  two  bi-trapezoidal 
plates  showed  that  the  developed  theory  agreed  quite 
closely  with  the  experimental  results.  The  discrep¬ 
ancies  between  experiment  and  theory  which  occurred  at 
large  longitudinal  curvatures  were  attributed  to  the 
initial  imperfections  in  the  plate. 
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-  one-half  the  plate  width 

-  thickness  of  the  longitudinal  edges 

-  modulus  of  elasticity 

-  normal  strains  in  the  x,y  and  z-directions 

-  one-half  maximum  plate  thickness  =  t  4-  c/2 

-  non-dimensional  coefficient  defined  in  figure  1.3 

-  bending  moments  per  unit  distance 

-  normal  forces  per  unit  distance 

-  shearing  forces  per  unit  distance 

-  radius  of  curvature  of  the  neutral  surface  in  the 
longitudinal  direction 

-  radii  of  curvature  of  the  middle  surface  in  the 
direction  of  the  subscript 

-  thickness  of  the  plate  at  a  point  x^ 

-  one-half  thickness  of  the  tapered  portion  of  the 
plate  as  shown  in  figure  1.3 


v  =  £*/  y 

w  -  deflection  of  the  midline  of  the  cross  section  in 

the  z-direction  measured  from  the  centroidal  axis 


w 

o 


initial  deflection  of  the  cross  section  measured 
from  the  centroidal  axis 


x3y5z  -  rectangular  coordinates  as  shown  in  figure  1 

a  =  2A(1+£)1/2 
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7 
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distortion  of  the  midline  of  the  cross  section 
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integral- 
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ent  solutions  of  the 


3(l-|i  ) 


1/2 


Rt 


O 
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INTRODUCTION 

When  a  beam  whose  width  is  the  same  order  of  mag¬ 
nitude  as  its  depth,  is  bent  into  a  uniform  longitudinal 
curvature  1/R,  the  transverse  curvature  is  -|a/R  producing 
an  anticlastic  surface.  This  is  provided  the  longitudinal 
edges  are  free.  When  the  width  of  the  beam  is  large  com¬ 
pared  to  its  depth  the  transverse  curvature  is  restrained 
in  the  central  portion  of  the  transverse  section,  however, 
the  edges  of  the  section  are  seen  to  curl.  The  bending 
behaviour  of  wide  beams  or  plates  has  received  consider¬ 
able  attention  in  the  literature.  Theoretical  work  has 
been  presented  by  Lamb^  and  Ashwell^2^  while  experiment¬ 
al  verification  of  Lamb's^  theory  has  been  given  by 
Bellow,  Ford  and  Kennedy ,  and  Conway  and  Nickola^  . 

These  investigators,  as  well  as  others,  have  shown  the 
mode  of  deformation  to  depend  on  the  ratio  b  /Rt  where  b 
is  the  beam  or  plate  width  and  t  is  the  thickness. 

The  problem  is  complicated,  however,  when  plates 
having  a  tapered  transverse  cross  section  are  considered. 
Plates  with  various  tapered  cross  sections  have  been  in¬ 
vestigated  theoretically  by  Flugge^  ,  Murray  and  Niles^, 

( 7 )  (8) 

Fung  and  Wittrickv  '  and  Conway  and  Farnham'  .  The  aim  of 
this  thesis  is  to  extend  the  work  of  the  above  authors  and 
to  consider  both  theoretically  and  experimentally  the  trans¬ 
verse  behaviour  of  plates  having  a  bi-trapezoidal  cross 


. 

. 
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section.  The  plate  orientation  and  a  general  bi-trapezoidal 
cross  section  are  shown  in  figures  1  and  2. 
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FIGURE  1.  PLATE  ORIENTATION 


FIGURE  2 .  GENERAL  BI-TRAPEZOIDAL  CROSS  SECTION 
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CHAPTER  I 

THEORETICAL  CONSIDERATIONS 
1.1  Historical  Review 

The  problem  of  finding  the  transverse  deflections 

of  a  plate  bent  by  uniformly  distributed  moments  on  its 

transverse  edges,  when  the  longitudinal  edges  are  free,  was 

(1) 

first  investigated  by  Lamb  in  1891.  His  development  of 
the  governing  differential  equation,  which  considered  the 
equilibrium  of  a  small  plate  element,  is  used  in  this 
thesis . 

In  analyzing  problems  with  similar  boundary  and 

(5) 

loading  conditions,  Flugge  m  1949,  derived  a  differ¬ 
ential  equation  which  accounted  for  a  thickness  variation 
in  the  transverse  direction.  This  work  included  results 
which  showed  the  transverse  deflections  for  rectangular, 
diamond  and  certain  bi-trapezoidal  sections  for  various 
values  of  longitudinal  curvature. 

(6) 

Also  in  1949,  Murray  and  Niles  developed  a 
graphical  method  and  its  numerical  equivalent  to  determine 
the  anticlastic  curvature  of  any  wide  beam  or  plate  of  an 
arbitrary  doubly  symmetric  solid  section.  This  numerical 
method  was  used  to  calculate  the  transverse  deflection 
curves  for  rectangular,  diamond  and  bi-convex  cross  sections. 

(7) 

Fung  and  Wittrick  in  1954  developed  the  govern- 
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ing  differential  equation  with  allowance  for  a  transverse 
thickness  variation  using  the  large  deflection  equations 
of  Von  Karman.  The  solution  of  the  differential  equation 

and  results  for  various  double-wedge  shapes  were  shown. 

(8) 

Recently  Conway  and  Farnhairr  '  ,  in  1964,  solved 
the  governing  equation  for  a  transverse  cross  section  which 
had  a  central  uniform  section  and  edge  sections  which 
tapered  in  a  parabolic  curve  to  form  an  infinitely  thin 
edge . 

The  governing  differential  equation,  which  is 
developed  below  for  the  transverse  deflections  of  a  bi- 
trapezoidal  plate,  is  completely  analogous  to  the  de¬ 
flections  of  a  cylindrical  tank  with  a  linearly  varying 
wall  thickness.  The  tank  problem  has  been  analyzed  by 
Timoshenko ^  and  Flugge^10^  . 

1.2  Development  of  the  Theoretical  Solution. 

1.2-1  Lamb's  Development 

Consider  the  plate  shown  in  figure  1.  It  is  sub¬ 
jected  to  a  longitudinal  curvature  produced  by  uniformly 
distributed  moments  applied  along  the  edges  which  lie 
in  the  x-z  plane.  The  x-z  plane  contains  the  transverse 
cross  section. 

The  following  assumptions  were  made: 


*  . 


\ 


. 
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1 .  The  behaviour  of  the  plate  under  loading  is  per¬ 
fectly  elastic  and  obeys  Hooke's  law, 

2.  The  material  is  isotropic  and  homogeneous. 

3.  The  slope  of  the  transverse  deflection  curve  is 
small  compared  to  unity. 

4.  The  weight  of  the  plate  is  negligible  compared 
to  the  applied  forces. 

Consider  the  equilibrium  of  the  plate  element  shown 
in  figure  1.1.  The  x  and  y-axes  are  chosen  to  be  principal 
axes  so  that  the  forces  and  moments  are  principal  valued. 

Force  equilibrium  parallel  to  the  x-axis  gives 

r  dN  ,  -] 

N  4-  ,_™2£  _  ^  dy  -  N  dy  =  0 
x  ox  2  x  2 

c>N 

therefore  ~  =  o  and  N  =  f  (y)  . 

ox  X  1 

At  the  edges  x  =  +b  there  are  no  surface  forces. 

It  follows  that  must  be  zero  at  these  edges  and  all 

across  the  transverse  section. 

N  =  0  1.2-1 

x 

Force  equilibrium  parallel  to  the  y-axis  shows  that 
must  be  a  constant  value. 


-  0 


dQ_ 


+ 


x 


bx 


dx 


dy 


Qxdy  + 


N  dxdO  =  0 

y 


f 
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ax 


FIGURE  1.1  PLATE  ELEMENT 


FIGURE  1.2  GEOMETRICAL  RELATIONSHIPS 
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Since  R^dQ  =  dy  this  reduces  to 


dQ 


x 


IX 


N 

+  H*  " 

y 


o. 


1.2-2 


Moment  equilibrium  about  the  y 


dM 


M  +  -v 
x  dx 


x 


dx 


axis  gives 
dQ 


dy  -  Mxdy  - 


Q  + 


x 


sir 


dx 


dxdy  =  0. 


Neglecting  the  second  order  differential  terms  as  being 
small  compared  to  the  first  order  the  result  is 


dM 


x 


Si  -  Q~  “  0 


1.2-3 


From  equations  1.2-2  and  1.2-3  it  follows  that 


d2M  N 

— =2  +  -X  -  0. 


1.2-4 


dx' 


R 


Hooke's  law  and  the  principle  of  superposition  applied  to 
a  plane  stress  problem  gives 


E 


(l-H2) 


(ex  +  (i€y)  and 


1.2-5 


E 


“y  M  2,  (ey+H£x)  • 

( 1~M-  1 

Considering  the  geometrical  relationships  shown  in  figure  1.2, 
the  total  strain  at  an  arbitrary  point  in  the  transverse 


cross  section  is 


OiJ  :•  'Of  e  .  J 


r:bxM  - 


prti scf  as  smi6^  I/cx^n.  .  sbso  bnoosa  s/i'j  pi 

leb^o  o;  axs<:  * 


2X  V3  if  :  ^rfi  i  :oq 
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I  I 

_  zZ  _  -  (z  -i-w)  _  -z  _  w 
ey  R  Ry  R-w  R-w 


Now 


M 


x 


t/2 


-t/2 


z  '  a  dz  1 
x 


1.2-6 


and  using  equations  1.2-5  and  1.2-6  this  becomes 

3 


M  = 
x 


-Et* 


1-  +  ±L- 

R  R 
x  y 


1.2-7 


12(l-|i  ) 

The  integration  for  the  moment  shows  that  the  first 
term  in  equation  1.2-6  is  due  to  bending.  The  second  term 
is  caused  by  stretching  or  contracting  the  middle  surface 
fibers . 

The  membrane  forces  induced  by  this  stretching  are: 

t/2  t/2 

N  =  /  o  dz '  and  N  =  /  a  dz '  . 

y  J  /  y  x  J  x 

-t/2  -t/2 


However  N  equals  the  stress  at  the  middle  surface  multi- 

X 

plied  by  the  thickness  t  therefore 


N  =  t 
x 


a 


x 


z  "  =0 


but  from  equation  1.2-1,  N  r  =  0 ,  which  gives 

X 


£x  + 


=  0. 


1.2-8 


z  '  -0 


r 
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Using  this  result 


N 

Y 


t/2 


a  dz  ' 

y 


t/2 


=  E 


-t/2 


rJz  '  = 


dz 


0 


-Ewt 

R-w 


1.2-9 


Substituting  equation  1.2-7  and  1.2-9  into  equation  1.2-4 
the  result  is 


Et' 


J_  + 

R  R-wl - 

x 


Ewt 
(R-w) 2 


0. 


Since  the  plate  is  bent  with  uniformly  distributed 
longitudinal  bending  moments,  each  transverse  section  is 
subjected  to  the  same  loading  conditions.  It  follows  that 
for  a  given  transverse  strip  w  is  a  function  of  x  only. 
Therefore,  the  partial  derivatives  in  the  above  equation 
can  be  replaced  by  total  derivatives. 

Also,  consider  R  where 
’  x 


1 


.2 

d  w 


dx' 


3/2 


-r—  <<C  1  then 
dx 

1  j;  d  w 

R  ,  2 
x  dx 


but  if 


'  •  l 
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The  resulting  equation  is 


Et' 


dx2  |l2(l-n2) 


d2w 


+ 


xir 


Ewt 


=  0. 


L  dx~  R-W  1  )  (R-w) 

Since  it  is  assumed  that  R-w  =  R  simplification  of  the 
above  gives 


dx 

(1) 


;3  dV)  +  U  (ty\  +  imzu!)  ^  .  o 

dx  y  R  dx  \  /  R^ 


1.2-10 


Lamb  solved  this  equation  by  considering  the  transverse 
section  to  be  of  constant  thickness. 


1.2-2  Extension  for  Transverse  Thickness  Variation 

(7) 

Fung  and  Wittrick  obtained  the  same  result  as 

( 2) 

Lamb  using  the  large  deflection  equations  of  von  Karman. 

(7) 

The  authors  considered  cases  where  the  thickness  varied 
along  the  cross  section  and  the  initial  shape  of  the  mid¬ 
line  was  distorted. 

Let  w  =  C  +  w  where  w  is  the  initial  distortion  of 
^  o  o 

the  midline  of  the  cross  section  measured  from  the  centroidal 
axis.  The  additional  deflection  due  to  the  applied  moments 
is  £.  Considering  wq  as  a  linear  function  of  x,  equation 
1.2-10  becomes 


dL  ft3  d iC]  +  U  d2(t3)  +  12(1-U2) 
dx2  V  dx2/  R  dx2  R2 


(£+wo) 


0. 


Consider  a  general  bi-trapezoidal  section  shown  in  figure  1.3 


snc  '  .  i  i.  >£  *eo 


IOX^  SX^bV  .  i;  'Y  )3'  .fiVS  I 


- 

.*  ex-xfiv  33  9m'x  ij  eriJ-  :i  fv.  e  a«o  b©’.  3ian  •• 

•  !,ir  -  ’J  t*  i,  -  .  ■  .*!  •  .  ;  -n  - 

-  •  • 

:  no  : ^  .tb 

o  in  Y  *i  T  -  it  iio*r  ss  it 
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FIGURE  1.3  TRANSVERSE  CROSS  SECTION 


The  thickness  t,  as  a  function  of  x,  is  given  by 

t  =  2tQ(l  -  - )  +  C  1.2-11 

for  any  positive  x.  If  the  edge  thickness  c  approaches 
zero  the  result  is  the  double-wedge  shape  considered  by 

(7) 

Fung  and  Wittrick  .  If  k  and  t  are  both  zero  the  cross 
section  becomes  rectangular. 

Equation  1.2-10  can  be  transformed  into  non-dim¬ 
ensional  terms  by  setting  p  =  t/t  and  £, ,  =  1  -  x/b  yield- 
ing 


1.2-12 


where 


. 
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By  making  the  additional  substitutions  (3  -  c/2tQ  and 
7  =  |  +  P,  the  thickness  ratio  is  given  by  p  =  27,  and 
equation  1.2-12  can  be  reduced  to 

(V  dftN  +  6^  +  a47(  =  0_ 

d7Z  \  dyV  R  0 

This  is  the  differential  equation  obtained  by  Fung  and 
( 7) 

Wittrick  for  double-wedge  cross  sections. 


1.2-3  General  Solution  of  the  Differential  Equation. 
Let  be  the  particular  integral  of  equation  1. 2-^13 

so  that 


4 

+  a  y£ 

so 


..4 

A  7w 


6nb2 

R 


7. 


This  gives 


w 

o 


1.2-14 


1.2-15 


The  total  solution  is 

5  “  C*  +  C0 

where  £  is  the  complimentary  function.  Equation  1.2-13 
then  becomes 


7 


3  ££ 
<37 2  . 


4  * 

+  aV  - 


0 


' 


' 
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and  division  by  7  results  in 


1  d2  /73  £]£\+  Ay  .  0_ 


7  d7‘ 


d7 


The  solution  of  this  equation  has  been  given  by  Timoshenko 
and  Flugge^*^  .  Equation  1.2-16  can  be  reduced  to  that  of 
two  equations  of  the  second  order  by  use  of  the  following 
identity. 


1.2-16 
(9) 


1  d“ 


7  d7‘ 


7 


3  dkl 

d72  - 


1  d_  ^  2  d_ 
7  d7  I  d7 


1  d  f  2 
~  ■ —  (  7 
L7  d7  \  d7  / 


T  .  .  .  1  d  f  2  df 

Letting  L(C  )  -  -  ^<7  ^ 


equation  1.2-16  reduces  to 


•k  j  4  * 

LU  )J  +  7s  X  =  0. 


1.2-17 


This  equation  can  be  written  in  one  of  the  following  forms. 

/  *  2  2  /  *  2 
L^L ( £  )  +  iAZ£  A  -  iA  )  +  iA  C  )  =  0 

L/L(C  )  -  i  A^£  )  +  iAZ  [L(^  )  -  iAZ^  J  =  0 


where  i  —  \/-l. 

Thus,  the  solutions  of  the  second  order  equations 


L(C  )  +  ix2c  =  0 


1.2-18 


must  be  solutions  of  1.2-17.  Since  these  equations  have 
an  imaginary  coefficient,  the  solutions  are  complex-valued 
functions  of  7,  and  the  solutions  of  one  equation  will  be 
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the  complex  conjugates  of  the  other. 

Since  the  solutions  are  complex  conjugates  it 
follows  that  they  are  linearly  independent  of  one  another. 
Hence,  the  two  independent  solutions  of  either  equation 
1.2-18  together  will  form  a  complete  system  of  four  in¬ 
dependent  solutions.  Assume 

q  “  ei  +  ie2 

?2  =  e3  +  i04 

are  two  independent  solutions  of 

*  p  * 

L(£  )  +  iAz£  =  0 

then 

h  '  ®1  -  iS2 

c4  =  03  -  ie4 

are  solutions  of 

L(£  )  -  iAZC  =  0. 

All  four  solutions  represent  the  complete  system  of  in¬ 
dependent  solutions.  By  using  sums  and  differences  of 
these  solutions  the  general  solution  of  1.2-17  can  be 
represented  as 

r  =  c^G.  +  c0©0  +  co0o  +  c.e. . 

D  Xi.  Z  Z 

To  determine  the  functions  0,  solution  of  only  one  of  the 


' 


•.  •'to  .  J  ;i  s 


rf:h  *(.i  »;.<  h 
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second  order  equations  above  is  needed.  Consider 

L(e  )  +  =  0. 

Writing  the  differential  operator  in  full  we  obtain 

y  +  2  ^  +  iA2C*  «  0. 

d72  d7 

Now  set  ip  *  2A(i7)1//2  =  rjs/iT  and  v  =  £v/y 
then  we  obtain 

, 2  d2v  ,  ,  dv  ,  /  ,  2  , v  _  « 
ip  r  f  +  if  -l)v  =  0. 

d^/  d^/ 

The  solutions  of  this  equation  are  the  Bessel  functions  of 
the  first  order  of  the  complex  argument  ip  giving 

v  =  AJ^ip)  +  B  H^)  . 

The  Bessel  functions  J ^ W >  H^(^)  of  the  first  order,  are 
connected  to  the  zero  order  functions  by 

Jl W  =  if  Jo(K'),  =  if  Ho(V-) 

The  real  and  imaginary  parts  of  Jo  and  Ho  may  be 
considered  as  real  functions  of  the  real  variable  r|.  These 
functions  are  known  as  Thomson  or  Kelvin  functions . 

The  Thomson  functions  are  introduced  by  the  fol¬ 
lowing  formulae 

Jo  (ip)  =  Jo(t]s/T)  =  berT)  -  i  bei  r| 

_  9 

Ho  (ip)  =  Ho(t]n/T)  =  — -  (kei  r\  +  i  ker  t}) 


. 


- 
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Differentiating  with  respect  to  f  and  separating  real  and 
imaginary  parts,  a  set  of  relations  for  the  first  order 
functions  are  obtained, 


(bei  '  r\ 


ber  '  t))  +  i(bei'r)  +  ber  '  iq) 


H.^) 


s/2 

TT 


(ker'r)  +  kei ' -q)  +  i  (ker'r) 


kei  '  r])  , 


where  the  prime  denotes  differentiation  with  respect  to  -q. 

The  solution  of  equation  1.2-16  can  be  written  as 


1 


S  ber'T)  +  T  bei  '  r\  +  U  ker '  iq  4-  V  kei'r) 


and  the  complete  solution  of  equation  1.2-13  becomes 


c 


1 

T1 


S  ber  1  t]  +  T  bei  '  r|  +  U  ker  '  r\  +  V  kei  '  t) 


-w 


6LLb< 

°  a4r 


1.2-19 


1.2-4  Constants  of  Integration 

At  the  free  edge  (x  =  b)  the  boundary  conditions 
show  that  the  moments  and  shear  forces  on  a  plane  whose 
normal  is  in  the  x  direction  are  zero  giving 


M  = 
x 


-Et" 


12(1  -  |±  ) 


d  w 

,  2 
dx 


+  1L_ 

Ry 


x=b 


V  = 
x 


Q  + 
x 


a 


^y  ^Mxy 


=  0 


'x=b 


0  and 


,  j  fr '  <9)1  V  H  ft' as*  U  +  ft'lri  1  +  _  •  vsd  8!  •  * 
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Since  M  =  0  the  condition  that  V  =0  reduces  to 
xy  x 


.3  -i 
d  w 

,  3 
dx  J 


=  0 


x=b 


Transformed  to  non-dimensional  terms  this  becomes 


d3;  3  a2;  ,  3  dc 

3  2  2 

■dr]  T)  drj  t)  dr)J 


r)=2AP 


1/2 


=  0. 


The  condition  =  0  at  the  free  edge  results  in 


,2 

d  w 


where 


L-dx  Jx®b  R 


R  -  w 


1.2-20 


and  transforming  this  to  non-dimensional  terms  gives 

u.b2  2” 

■V2  "  '  4A4H  ^  — 1/2  ' 


rd2C  1  d£ 
2 


Ldr)  ri  dr|-lTf::2Af3 


rj2  A|3' 


The  symmetry  of  the  problem  allows  the  condition 
In  non-dimensional  terms  this  is 


dx 


1.2-21 


-  0. 


x=0 


K 

_dTL 


n==2  A  ( l+f3) 


1/2 


=  0. 


1.2-22 


The  fourth  condition  results  from  the  boundary  condition 

b 
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and  can  be  simplified  with  equation  1.2-9  giving 

b 

f  (w  +  P)  tdx  =  0 . 

J  ° 

o 

Since  the  problem  is  symmetrical  about  the  longit¬ 
udinal  axis,  this  integral  is  written  for  only  half  the 
transverse  cross  section.  The  transformation  of  the  above 
equation  to  non-dimensional  terms  gives 

1+P 

£ydy  =  0.  1.2-23 


Before  evaluating  equations  1.2-20  to  1.2-23  wq  must  be 
determined.  For  the  section  considered  in  figure  1.3,  wq  is 
given  by 


w 


=  -t  (l+kW  +  t  (1+k)  +2/3J 


<2(3  +  1) 


The  four  conditions  given  by  equations  1.2-20  to  1.2-23 
result  in  four  equations  for  the  constants  S,T,U  and  V. 
These  are  given  below  where,  for  simplification, 


a  =  2A(l-tf3) 


1/2 


2Ap 


1/2 


0  =  S 


-e"3  ber  e  -  24e  bei  e  -  48  ber '  e  +  8e^  bei'e 


+  T 


bei  e  +  24e  ber  e  -  48  bei'e  -  8e^  ber 'e 


+  U 


+  V 


■e^  her  e  -  24e  kei  e  -  48  ker '  e  +  8e^’  kei'e 


-e^  kei  e  +  24e  ker  e  -  48  kei'e  -  8e^  ker ' e 


1 . 2- 20a 


' 


20 


0  =  S 

+  T 

+  U 

+  V 


4e  bei  e  -  e  bei'e  +  8  ber '  e 


-4e  ber  e  +  e  ber ' e  +  8  ber ' e 


4e  kei  e  -  e  kei'e  +  8  ker '  e 


■4e  ker  e  +  e  ker ' e  +  8  kei'e 


+ 

RA2 


1 . 2-21a 


0  =  s 

-abei  a  -  2  ber 'a 

+  T 

a  ber  a  -  2  bei 'a 

+  u 

-akei  a  -  2  ker 'a 

—  — 

+  V 

a  ker  a  -  2  kei ' a 

■f 


t.  a  ( l  fk) 
__o _ _ 

2  A^ 


1 . 2-22a 


0  =  S 


2  2 
a  ber  a  -  2a  bei' a  -  e  ber  e  +  2e  bei'e 


+  T 


2  2 
a  bei  a  +  2a  ber 'a  -  e  bei  e  -  2e  ber ' e 


+  U 


a 


2  2 

ker  a  -  2akei'  a  -  e  ker  e  +  2e  kei'e 


+  V 


2  2 
a  kei  a  +  2aker '  a  -  e  kei  e  -  2e  ker ' e 


6ub‘ 


a4r 


,4  4V 
~  e  ) 


1.2-23a 


The  solution  of  the  above  equations  for  the  constants 
in  the  general  solution  was  obtained  using  an  IBM  7040/1401 
digital  computer.  Details  of  this  solution  appear  in  Appendix  I. 


- 


,  .  -  l  ••  •  '  -  -  •  ■'  -  -  ■  ! 
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1.3  Comparison  of  Results  to  Previous  Work. 

Earlier  investigators  showed  solutions  for  a  uni¬ 
form  plate  and  for  a  tapered  plate  with  a  line  edge. 

(12) 

Bellow'  gave  the  results  of  Lamb's  solution  for  a  uni- 

( 7) 

form  plate.  Fung  and  Wittrick  gave  results  for  the 
linearly  tapered  plate  with  a  line  edge.  By  changing  para¬ 
meters  in  the  above  solution,  both  of  these  earlier  results 
can  be  closely  approximated. 

Using  the  theory  developed  above,  neither  the  flat 
plate  nor  the  tapered  plate  with  line  edge  solution  can  be 
reached  exactly.  For  the  uniform  plate,  the  thickness  of 
the  tapered  portion,  2tQ,  must  be  zero.  This  would  cause 
the  value  of  A  to  become  infinite.  The  tapered  plate  with 
the  free  edges  being  lines  is  obtained  when  the  edge  thick¬ 
ness  c  goes  to  zero.  With  c  zero,  the  variable  tj  is  zero 
at  the  free  edges  and  the  ker '  function  becomes  negative 
infinity.  For  this  case  of  zero  edge  thickness,  the 
boundary  conditions  of  zero  moment  and  zero  shear  along 
the  free  edges  no  longer  have  any  meaning.  Even  though 
these  solutions  can  not  be  reached  exactly  they  can  be 
closely  approximated  by  giving  tQ  and  c  values  which  are 
small  compared  to  the  total  thickness  of  the  plate. 

In  figure  1.4  the  deflections  of  a  uniform  plate 
are  compared  to  those  using  the  theory  developed  above.  These 
curves  were  obtained  for  a  value  of  t  which  was  approxi- 


' 


SOLUTION  FOR  b  /Rd  =  10. 


22 


S-J1 


I 


I 


I 


FIGURE  1.4  COMPARISON  OF  THEORETICAL  RESULT  WITH  LAMB’S  SOLUTION 
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X\A 


FIGURE  1.5  COMPARISON  OF  THEORETICAL  RESULT  WITH  FUNG  AND  WITTRICK'S  SOLUTION 
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mately  4  percent  of  the  total  thickness.  The  curves  are 

2 

compared  on  the  basis  of  the  same  b  /Rt  ratio.  Actual 

(12) 

flat  plate  results  were  obtained  from  work  done  by  Bellow 

2 

since  these  were  available  for  a  large  range  of  b  /Rt  ratios. 
It  must  be  noted  that  the  b  used  by  Bellow  was  the  total 
plate  width,  not  half  the  width  as  used  in  the  theoretical 
development  of  section  1.2. 

Figure  1.5  shows  the  comparison  to  Fung  and  Wittrick' 
results.  Their  results  were  taken  from  graphs  so  the  re¬ 
production  in  figure  1.5  may  not  be  exact.  To  compare  with 
these  results  a  solution  for  c  equal  to  approximately  4  per¬ 
cent  of  the  maximum  thickness  is  used.  The  definition  of  A 

(7) 

is  the  same  as  used  by  these  authors 

The  agreement  obtained  in  each  of  the  above  cases 
is  good.  In  the  following  chapters,  experimental  verif¬ 
ication,  for  specific  cases  which  lie  between  the  extremities 
discussed  above,  is  undertaken. 


(7) 
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CHAPTER  II 


EXPERIMENTAL  CONSIDERATIONS 


2„1  Numerical  Integration  Technique 

The  curvature  in  the  transverse  direction  due  to  the 
distortion  of  the  midline  is 


2 

1  _  d  ^ 

R  .2 
x  dx 


when  using  the  assumption  that  the  square  of  the  slope, 


is  small  compared  to  unity, 
this  expression  with  respect  to  x, 


transverse  direction,  results  in 


Double  integration  of 
the  coordinate  in  the 


C(x)  = 


dxdx 


where  c^  and  c^  are  arbitrary  constants. 

Consider  a  transverse  strip  with  curvature  1/R 

i 

at  the  point  x^.  Let  the  strains  indicated  by  the  two 

strain  gauges  mounted  on  the  upper  and  lower  surfaces  of 

this  plate  at  the  point  x^  be  e^  and  e^  ,  respectively. 

i  i 

The  bending  strain  at  a  point  in  a  fiber  is  known  to  be 
directly  proportional  to  the  curvature  at  that  point  by 
the  relation 

z 


e  = 


where  z  is  the  distance  from  the  neutral  surface  to  the 
point.  Considering  the  upper  surface  of  the  plate,  the 
bending  strain  becomes 
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^ebending) 


t . 
 1 


2R 


2.2-1 


x 


where  t^  is  the  thickness  of  the  strip  at  x^. 

Since  the  actual  measured  strains  may  include  mem¬ 
brane  strains  these  must  be  subtracted.  The  membrane  strain 


is  given  by  the  average  value  of  ern  and  eT,  with  the  re- 

1  .  B  ♦ 

1  1 

suiting  bending  strain 


(ebendingt  eT. 


(eT.  +  eB.> 

1  1 


eT.  “  eB. 
1  1 


and  the  curvature  is  then  given  by 


2.2-2 


l 

To  calculate  the  deflections  from  the  curvatures 

required  a  double  integration  which  was  done  numerically. 

(13) 

The  trapezoidal  rule  which  was  used  m  this  process 

(3) 

has  been  described  by  Bellow  et  al  .  The  first  integ¬ 
ration  results  in  the  slope  of  the  distortion  curve.  This 
distortion  curve  is  symmetrical  about  the  center  of  the 
transverse  section  so  the  integration  was  started  at  the 
point  x  =  0  where 


•  r  r  k 
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dx 


x=0 


0. 


The  second  integration,  which  gave  the  distortion  of  the 
midline  of  the  section,  was  also  started  at  the  plate 
center.  The  distortion  of  the  midline,  was  calculated 
relative  to  the  center  of  the  cross  section  (x  =  0) .  To 
comply  with  the  boundary  conditions  in  Chapter  I  the  dis¬ 
tortion  was  to  be  measured  relative  to  the  centroidal  axis, 
so  the  trapezoidal  rule  was  used  again  to  calculate  the 
position  of  the  centroid. 

If  the  distance  between  two  neighbouring  measuring 
stations  is  h  =  -  x^  then  the  length  of  the  line  join¬ 

ing  the  points  representing  the  distortion  at  x^  and  x^+1 


+ 


The  first  moment  of  area  of  this  portion  of  the  cross  sect¬ 
ion  (between  x^  and  x^_^)  about  the  horizontal  axis  through 


=  0 


x=Q 


is 


+^i+i 


2.2-3 


Multiplying  equation  2.2-3  by  h/h  the  result  is 


h 


<■  +  C- 


i+1 


+  fci+r 


1  + 


' 
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Assuming  again  that  the  square  of  slope  of  the  dis¬ 
tortion  curve  is  small  compared  to  unity,  the  last  term  in 
the  above  is  negligible.  The  moment  of  this  segment  of 
the  cross  section  reduces  to 


The  centroid  of  the  distorted  curve  is  given  by 

b 

Za  =  I"  rtdx 

-b 

where  A  is  the  area  of  the  cross  section.  Therefore,  £  is 


(«o  +  +  0 


+  t 


n 


/ 


2.2-4 


where  n  is  the  number  of  measuring  stations.  The  more 
measuring  stations  used  the  closer  the  above  analysis  ap¬ 
proaches  the  actual  case.  The  distortion  of  the  midline  is 
given  by  CTrue  =  C  ~  Z  and  the  position  of  the  midline  of 
the  cross  section  is  given  by 


w  =  £ 


m  +  w 
True  o 


2.2-5 


The  initial  deflection  of  the  midline  of  the  cross 


section,  w  ,  is  given  relative  to  the  centroidal  axis. 
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2,2  Experimental  Apparatus 

2,2-1  Plates 

The  plates  tested  were  machined  from  a  flat  rolled 

•v 

ie 

sheet  of  Alclad  7Q75-T6  .  This  material  has  a  modulus  of 
elasticity  of  10,4  x  10^  psi  and  a  Poisson's  ratio  of  0.333 
and  was  chosen  because  it  had  a  higher  yield  stress  and  a 
lower  modulus  of  elasticity  than  mild  steel.  The  high 
yield  stress  allowed  larger  longitudinal  curvatures  than 
possible  with  mild  steel.  With  a  lower  modulus  of  elast¬ 
icity  than  steel  a  specific  curvature  could  be  obtained 
with  a  smaller  load. 

The  first  plate  tested,  Plate  I,  was  machined  from 
a  sheet  of  aluminum  0.125  x  12.0  x  67.0  inches.  The  bi- 
trapezoidal  cross  section  was  obtained  using  an  end-mill 
cutter  on  a  vertical  milling  machine.  Figure  2.1  shows 
half  the  transverse  cross  section  of  the  finished  plate 
and  the  cross  section  as  it  was  assumed  to  be  for  calcul¬ 
ation  purposes.  During  the  machining  a  warping  of  the 
entire  cross  section  occurred  (the  bottom  edge  was  initially 
flat) ,  with  pronounced  curling  at  the  edges.  The  cross 


*  Designation  and  properties  of  the  alloy  used  were  ob¬ 
tained  from  Alcoa  Structural  Handbook,  Aluminum  Company 
of  American,  1956,  p35. 
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FIGURE  2.1  COMPARISON  OF  ACTUAL  AND  ASSUMED  CROSS  SECTION  FOR  PLATE 
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FIGURE  2 . 2  o  COMPARISON  OF  ACTUAL  AND  ASSUMED  CROSS  SECTION  FOR  PLATE  II 


section  dimensions  were  chosen  so  that  A  could  be  varied 
over  the  range  0<-A<i7.5. 

Plate  II  was  cut  from  aluminum  sheet  0.25  x  8.0  x  70.0 
inches.  The  final  section,  which  is  shown  in  figure  2.2 
was  obtained  using  a  horizontal  milling  machine  with  a  4- 
inch  long,  2.5-inch  diameter  cutting  tool.  This  tool 
allowed  one  side  of  the  section  to  be  machined  without  moving 
the  tool  transverse  to  the  plate.  The  result,  shown  in 
figure  2.2  shows  the  straight  edges  obtained  by  this  method. 
The  section  warped  during  the  machining  but  this  warping 
was  quite  linear,  allowing  an  accurate  description  of  the 
actual  section.  As  a  result  of  the  final  shape  of  the  cross 
section,  A  values  up  to  approximately  4.0  were  all  that 
could  be  obtained  experimentally  without  yielding  the  mat¬ 
erial  . 


2.2-2  Strain  Gauges 

Electrical  resistance  strain  gauges  were  mounted 
along  the  midline  on  the  upper  and  lower  surfaces  of  the 
plate  in  the  transverse  direction.  The  position  of  the 
gauges  on  Plates  I  and  II  are  shown  in  figures  2.3  and  2.4, 
respectively,  where  the  bottom  gauges  were  mounted  directly 
below  those  on  top.  In  order  to  keep  the  error  in  using 
the  trapezoidal  rule  as  small  as  possible,  a  large  number 
of  measuring  stations  should  be  used.  A  discussion  of  the 
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FIGURE  2.4  PLATE  II 
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truncation  error  is  given  in  Chapter  III.  It  was  decided 
that  13  stations  on  half  the  cross  section  would  give  the 
required  accuracy  while  the  application  of  this  number  of 
gauges  would  not  be  difficult.  As  a  result  the  distance 
between  gauge  centers  was  0.5  inches  for  Plate  I  and  0.333 
inches  for  Plate  II.  To  obtain  the  gauge  center  distance 
of  0.333  inches  for  Plate  II  the  gauges  were  staggered. 

As  well  as  these  transverse  gauges,  four  longitudinal  gauges 
on  each  plate  (two  on  top  and  two  opposite  on  the  bottom) 
were  used  to  calculate  the  longitudinal  curvature. 

The  strain  gauges  used  were  “Budd"  metal  film 
electrical  resistance  gauges  type  C12-121.  These  gauges, 
which  are  temperature  compensated  for  aluminum,  have  a  grid 
length  of  0.125  inches  and  an  overall  length  of  0.250  inches. 
Two  advantages  of  these  gauges  are  the  thin  backing  material 
on  the  gauge  and  low  transverse  sensitivity.  The  thinner 
the  backing  material  the  closer  the  actual  gauge  grid  is 
to  the  surface  of  the  plate.  This  allows  a  closer  measure¬ 
ment  of  the  strain  at  the  surface.  Low  transverse  sen¬ 
sitivity  is  important  since  the  plate  was  under  large 
strains  in  the  longitudinal  direction  compared  to  those 
measured  in  the  transverse  direction. 

To  obtain  the  strains  as  close  to  the  edge  of  the 
plate  as  possible,  type  C12-111  gauges  were  used.  These 
have  a  grid  length  of  0.063  inches  with  the  result  that  the 
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gauge  center  was  0.0315  inches  from  the  edge.  Since  the 
distance  between  gauges  was  assumed  to  be  constant  the 
small  error  introduced  was  assumed  to  be  negligible  in 
the  analysis. 


2.2-3  Loading  Apparatus 

The  method  of  applying  a  uniformly  distributed  mom 

ent  to  the  longitudinal  ends  of  the  plate  was  similar  to 

(4) 

that  used  by  Conway  and  Nickola  .  The  two  ends  of  the 
test  plate  were  bolted  on  one-inch  diameter  shafts.  These 
shafts  were  keyed  to  pulleys,  and  loading  of  the  plates  was 
done  by  applying  loads  to  the  rims  of  these  pulleys.  This 
loading  system,  which  was  supported  on  bearings  to  reduce 
the  frictional  effects,  is  shown  in  figure  2.5. 

The  loading  arrangements  discussed  above  had  the 
advantage  of  being  adaptable  to  different  length  and  width 
of  plates.  Pulley  sizes  were  easily  changed  allowing 
reasonable  loads  to  be  used  in  obtaining  the  desired  curv¬ 
ature.  A  6-inch  diameter  pulley  was  used  with  Plate  I 
while  a  15-inch  diameter  pulley  was  used  with  Plate  II. 

To  check  the  influence  of  the  support  condition  on 
the  tapered  plate,  a  second  set  of  transverse  strain  gauges 
were  mounted  on  Plate  I  approximately  one  foot  from  one 
of  the  edges  being  loaded.  The  results  of  tests  with 
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FIGURE  2.5  VIEW  OF  APPARATUS 
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these  gauges  showed  little  or  no  difference  to  those  loc¬ 
ated  2.5  feet  from  the  loaded  edges.  Other  tests  were 
made  with  the  supports  each  18  inches  from  the  section 
being  tested.  Again,  little  difference  in  the  results  was 
observed. 


2.2-4  Recording  and  Associated  Equipment 

■it 

A  multichannel  digital  data  processor  was  used  to 
measure  the  voltage  drop  across  the  strain  gauges  and  con¬ 
vert  this  analog  reading  to  a  digital  one.  This  apparatus 
was  calibrated  so  that  the  output  was  the  strain  measured 
in  microinches  per  inch.  Recording  of  the  strains  was 
done  with  an  IBM  024  card  punch.  Since  a  card  punch  had 
not  been  used  in  conjunction  with  the  processor  previously, 
a  description  of  the  electrical  hook-up  of  the  card  punch 
to  processor  is  given  in  Appendix  III.  A  picture  of  this 
system  is  shown  in  figure  2.5. 

The  strains,  which  were  recorded  on  IBM  cards,  were 
used  as  data  for  a  Fortran  IV  computer  program.  This 
program,  which  was  run  on  the  IBM  7040/1401  system,  cal- 
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Design  and  operation  of  this  processor  has  been  discussed 

by  Bellow  (12) . 
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culated  the  deflections  by  the  numerical  method  outlined 
in  section  2.2.  Appendix  II  includes  a  description  of  the 
program  used. 


2.3  Test  Procedure 

The  procedure  for  testing  both  plates  was  essentially 
the  same.  The  plate  being  tested  was  initially  put  on  a  flat 
surface  where  the  "zero"  or  initial  readings  of  the  strains 
were  taken.  Actually  the  zero  balancing  system  of  the 
processor  was  used  to  bring  the  strains  as  close  to  zero  as 
possible.  The  "zero"  readings  were  then  taken.  These  initial 
strains  were  subtracted  from  those  taken  under  load  to  ob¬ 
tain  the  absolute  change  in  strain  due  to  the  load.  The 
plate  was  loaded  by  adding  weights  to  the  loading  hangers 
and  thereby  applying  a  uniformly  distributed  moment  to  the 
transverse  edges  of  the  plate.  At  the  various  A  values  the 
transverse  strains  were  recorded  on  the  card  punch.  The 
output  of  the  processor  was  also  printed  on  a  typewriter  to 
give  a  check  on  the  strain  readings.  If  the  readings  were 
abnormal  they  could  be  disregarded  and  another  set  taken. 

Once  the  largest  curvature  was  reached,  and  the  transverse 
strains  recorded,  the  plate  was  unloaded.  During  the  un¬ 
loading  the  transverse  strains  were  again  recorded  at  the 
various  A  values.  After  unloading,  with  the  plate  again 
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resting  on  a  flat  surface,  another  set  of  strains  were  re¬ 
corded  and  these  checked  against  the  initial  readings  to 
see  if  the  strain  gauges  were  functioning  properly. 

It  was  found,  that  after  a  loading  program,  as  in¬ 
dicated  in  the  above  paragraph,  some  of  the  gauges  tended 
to  drift.  In  order  to  minimize  this  drift,  it  was  decided 
to  take  a  set  of  "zero"  strain  readings  before  each  load 
was  applied  to  the  plate.  These  "zero"  strains  were  sub¬ 
tracted  from  the  strains  recorded  under  the  load  next 
applied  to  give  an  absolute  strain  reading.  More  consis¬ 
tent  results  were  obtained  by  this  testing  scheme  than  by 
the  first  procedure. 

As  many  as  10  tests  were  taken  at  each  A  value.  The 
results  of  these  tests  are  shown  in  the  following  chapter. 
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CHAPTER  III 


EXPERIMENTAL  AND  THEORETICAL  RESULTS 


3.1  Remarks  on  Experimental  Results 

Since  the  plates  and  the  loads  applied  to  them  were 
symmetrical  about  the  longitudinal  y-axis,  the  results  of 
both  halves  of  the  transverse  cross  section  were  compared 
to  the  theoretical  curves  for  half  the  cross  section.  In 
this  way  two  sets  of  experimental  results  were  obtained 
from  each  test.  The  plates  were  loaded  to  produce  the 
longitudinal  curvatures  corresponding  to  the  various  A 
values  at  least  5  and  sometimes  10  times  each.  As  a  re¬ 
sult,  there  were  at  least  10  and  up  to  20  sets  of  results 
to  compare  for  each  A  value. 

Figure  3.1  shows  a  typical  range  of  values  ob¬ 
tained  experimentally  for  Plate  I.  When  scattering  of 
results  occurred  as  shown  in  figure  3.1,  more  than  the 
minimum  of  five  tests  were  taken.  The  results  shown  in 
the  following  sections  for  Plate  I,  are  the  averages 
obtained  from  a  series  of  tests  at  the  particular  A  value. 

The  experimental  results  of  Plate  II  were  more 
precise  than  those  of  Plate  I.  The  graphical  representation 
of  these  tests  showed  the  points  from  various  tests  to  be 
almost  indistinguishable  from  one  another.  The  results  of 
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between  five  and  seven  tests  were  averaged  and  compared  with 
the  theoretical  curve  for  a  specific  A  value. 

3.2  Comparison  of  Experiment  and  Theory 

Initially  it  was  felt  that  the  warping  of  the  plates , 
which  occurred  during  the  machining,  would  not  affect  the 
final  results.  After  testing  Plate  I  extensively  it  was 
found  that  the  agreement  between  theoretical  and  experi¬ 
mental  results  was  poor  if  the  lower  edge  of  the  section 
was  assumed  to  be  flat  (k  =  0) .  It  was  decided  to  check 
the  bending  stresses  in  the  transverse  direction  to  see 
if  there  were  any  errors  in  the  numerical  integration 
technique  used  to  obtain  the  experimental  deflections. 

The  bending  stress  in  the  x-direction  at  point 
x.  of  the  transverse  cross  section  is  given  by 


3.2-1 


are  the  bending  strains  in  the  x 


1 


1 


and  y  directions,  respectively.  These  bending  strains  are 


given  by  equation 


and 
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x . 
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where  again  t ^  is  the  plate  thickness  at  x.. . 


Using  the 
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assumption  that 


1_  i  ill 

R  j  2  ’ 

x  dx 


then 
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and  since 


R  -  w  = 
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3.2-3 


Substituting  equations  3.2-2  and  3.2-3  into  3.2-1  the 
result  is 
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This  formula  was  used  to  calculate  the  transverse 
bending  stress  theoretically.  The  experimental  bending 
stresses  were  also  found  using  equation  3.2-1.  The  value 


of  (  e 


'y 


is  the  same  as  used  in  the  theoretical  result 


B'  . 
l 


in  the  above  paragraph.  The  bending  strain  in  the  x- 
direction  ( e  ;  was  found  to  be 
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in  section  2.2.  The  bending  stress  (  o  from  the  ex 
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periment  was  then  calculated  from 
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B/i  2(l-n) 
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3.2-6 


where  R  was  measured  by  the  longitudinal  strain  gauges. 


The  results  of  these  "stress  checks"  showed  that 


the  experiments  and  theory  still  deviated  considerably 
and  that  the  numerical  integration  was  not  the  entire 
reason  for  the  discrepancies  in  the  deflection  curves. 

Examination  of  equation  3.2-5  showed  the  only 
variable  for  a  specific  plate  was  k,  a  measure  of  the 
amount  of  warping  of  the  section.  Using  values  of  k  in 
the  range  0z~k<1.0  agreement  between  theory  and  experi¬ 
ment  for  the  bending  stresses  was  obtained.  The  trans¬ 
verse  cross  section  of  Plate  I  being  tested  was  measured 
and  it  was  found  that  the  warping  at  the  center  line 
corresponded  to  a  value  of  k  =  0.909.  The  bending 
stresses  obtained  from  the  experiments  were  compared  to 
theoretical  stresses  calculated  for  k  ~  0.0  and  k  =  0.909. 
The  two  results  gave  closer  agreement  assuming  k  =  0.909. 
In  a  similar  manner  Plate  II  was  measured  and  the  warp¬ 
ing  found  to  correspond  to  a  k  =  0.727.  A  check  of  the 
bending  stresses  showed  that  the  theoretical  result  for 
k  =  0.727  gave  closer  agreement  to  the  experiments  than 
k  =  0.0.  Figure  3.2  and  3.3.  show  typical  comparisons 


between  the  theoretical  and  experimental  bending  stresses. 
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FIGURE  3.2  STRESS  CHECK  FOR  PLATE  I  AT 
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As  a  result  of  these  "stress  checks",  the  trans¬ 
verse  deflection  curves  shown  in  figures  3.4  to  3.11  for 
Plate  I  and  in  3.12  to  3.15  for  Plate  II  use  k  =  0.909  and 
k  =  0.727,  respectively. 

3.3  Discussion  of  Results 

Figures  3.4  to  3.15  indicate  that  good  agreement 
between  experiment  and  theory  was  obtained  for  Plate  II 
and  for  the  lower  A  values  of  Plate  I.  The  disagreement 
seen  at  A>B .  0  for  Plate  I  could  have  been  caused  by  some 
of  the  following  sources  of  errors 

(1) ,  error  in  calibrating  the  digital  data  processor 

(2)  error  in  measurement  of  the  strains 

(3)  error  in  using  the  trapezoidal  rule. 

(12) 

Bellow'  has  shown  that  the  calibration  error  for 
the  processor  was  approximately  + 0 . 8  percent,  and  the 
measurement  error  was  +0.6  percent  full  scale.  Since  at 
the  higher  A  values  the  strains  recorded  were  in  the 
order  of  one-quarter  to  one-half  full  scale,  the  measure¬ 
ment  error  was  not  more  than  approximately  +2.5  percent. 

The  truncation  error,  which  resulted  from  using  the 

(12) 

trapezoidal  rule,  has  been  shown  by  Bellow  to  be  m 

the  order  of  +0.25  percent  when  using  a  one-eighth  inch 
plate  with  measuring  stations  one-half  inch  apart.  The 
total  error  from  the  above  causes  was  less  than  4.0  per- 
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cent  and  applying  this  error  to  the  values  obtained  made 
little  change  in  the  experimental  results.  The  dis¬ 
crepancies  in  the  curves  for  Plate  I  for  A>3.0  were  due 
to  much  larger  errors  than  those  discussed  in  this  para¬ 
graph  . 

It  was  believed  that  the  reason  the  deflections 
of  Plate  I  obtained  experimentally  differed  from  the 
theoretical  result  was  the  non-uniform  warping  which 
occurred  during  machining.  Examination  of  figure  2.1 
shows  that  the  assumed  and  actual  cross  sections  of 
Plate  I  deviated  considerably,  while  figure  2.2  shows 
that  the  assumed  cross  sectional  shape  of  Plate  II  closely 
approximated  the  actual  case.  Since  good  agreement  was 
obtained  from  tests  on  Plate  II  it  was  concluded  that  the 
discrepancies  found  in  Plate  I  were  due  to  non-uniform 
warping  of  the  transverse  cross  section. 


3.4  Conclusions 

l.The  theory  presented  in  Chapter  I  predicted  with 
reasonable  accuracy,  the  experimental  results  of 
the  two  plates  tested  in  the  range  0-=^A<r4.  The 
discrepancies  seen  at  the  higher  A  values  for 
Plate  I  were  probably  due  to  the  initial  imper¬ 
fections  of  this  plate. 
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2.  The  earlier  theoretical  results  of  Lamb^^  for 

a  rectangular  transverse  cross  section,  and  Fung 

(7) 

and  Wittnck  for  a  double-wedge  cross  section 
were  shown  to  be  limiting  cases  of  the  theory  of 
Chapter  I . 
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FIGURE  3.6  TRANSVERSE  DEFLECTIONS  FOR  PLATE  I  AT  A  =  2 . 5 


1.0  J 

FIGURE  3.7  TRANSVERSE  DEFLECTIONS  FOR  PLATE  I  AT  A  -  3 . 0 
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FIGURE  3.8  TRANSVERSE  DEFLECTIONS  FOR  PLATE  I  AT  A  = 
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FIGURE  3-11  TRANSVERSE  DEFLECTIONS  FOR  PLATE  I  AT  A  ~  7.5 
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FIGURE  3-12  TRANSVERSE  DEFLECTIONS  FOR  PLATE  II  AT  A  =  1.0 
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FIGURE  3-13  TRANSVERSE  DEFLECTIONS  FOR  PLATE  II  AT  A 
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FIGURE  3-14  TRANSVERSE  DEFLECTIONS  FOR  PLATE  II  AT  A  = 
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FIGURE  3-15  TRANSVERSE  DEFLECTIONS  FOR  PLATE  II  AT  A  =  3.8 
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CHAPTER  IV 

CONCLUDING  REMARKS 

4.1  Summary  of  the  Thesis 

In  Chapter  I  theory  governing  the  transverse  de¬ 
flections  of  thin  plates  with  a  bi-trapezoidal  cross 
section  was  given.  This  theory  was  developed  by  extend¬ 
ing  the  work  of  Lamb^  and  Fung  and  Wittrick^^  .  The 

developed  theory,  which  lies  between  the  two  extremes  dis- 

(1  7) 

cussed  by  the  authors  9  ,  was  shown  to  approximate  these 

earlier  results  by  a  suitable  selection  of  the  variable 
parameters . 

With  the  aid  of  a  digital  computer,  the  theory 
derived  in  Chapter  I  was  used  to  predict  the  experimental 
results  of  tests  carried  out  on  two  bi-trapezoidal  plates. 
The  computer  was  again  used  to  integrate  the  strains, 
measured  on  the  upper  and  lower  surfaces  of  the  plates,  to 
obtain  the  deflection  curves  of  the  transverse  cross  sect¬ 
ional  midline. 

The  agreement  between  the  predicted  and  actual  re¬ 
sults  was  good  for  Plate  II  and  for  Plate  I  at  the  lower 
A  values.  The  disagreement  at  values  A>3.0  for  Plate  I 
was  thought  to  be  caused  by  the  non-linear  warping  of  the 
cross  section.  The  discrepancies  between  experiment  and 
theory  are  most  evident  at  the  edges  of  the  cross  section 
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where  pronounced  distortion  occurred  during  machining. 


4.2  Practical  Applications  and  Problems  for  Further  Study 

The  theory  developed  above  is,  of  course,  directly 
applicable  to  other  bi-trapezoidal  plates  with  similar  bound¬ 
ary  conditions.  This  type  of  problem  occurs  in  thin  wing 

(5) 

design  as  pointed  out  by  Flugge  .  In  this  case  the 

transverse  deflections  can  influence  the  bending  stiffness 

as  well  as  the  aerodynamic  qualities  of  the  wing.  Conway 

and  Farnhamv  '  have  discussed  the  transverse  distortion 

problem  in  connection  with  magnetic  tapes  used  in  computer 

applications.  The  tapes  when  bent  over  a  reading  head 

distort  causing  uneven  wear.  Using  a  tape  cross  section 

other  than  rectangular  could  reduce  this  distortion. 

The  theoretical  moment-curvature  relationship  for 

(5) 

certain  bi-trapezoidal  plates  has  been  shown  by  Flugge 

This  could  be  extended  to  all  bi-trapezoidal  plates  and  the 

instability  in  the  applied  moment  similar  to  that  shown 

( 14) 

theoretically  by  Ashwell  for  uniform  thickness  plates, 

could  be  investigated  theoretically  and  experimentally. 
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APPENDIX  I 


FORTRAN  IV  PROGRAM  FOR  EVALUATION  OF  DEFLECTIONS  AND 

BENDING  STRESSES . 

In  Chapter  I  the  solution  of  the  differential 
equation  for  the  transverse  distortion  of  a  bi-trapez- 
oidal  plate  was  shown  to  be 


S  ber '  T)  +  T  bei  '  rj  +  U  ker  '  rj  +  V  kei'r) 


-  w 


.  Sub' 


4 
A  R 


The  program  shown  below  solves  equations  1.2- 20a 
to  1.2-23a  for  the  constants  S,T,U  and  V.  This  is  accomp¬ 
lished  through  the  use  of  a  subroutine  called  SUBROUTINE 
SOLVE  (shown  below) .  These  values  of  S,  T,  U  and  V, 
called  X(l)  to  X(4)  in  the  program  are  used  in  calculating 
the  distortion  The  deflection  of  the  midline  of  the 

cross  section,  w  =  T  +  w  ,  is  calculated  and  made  non- 

7  Q  7 

dimensional  by  dividing  by  h  half  the  total  maximum 
thickness  of  the  plate.  Finally  the  bending  stresses, 
called  SX( J)  are  computed  using  equation  3.1-5. 

The  values  of  the  Thomson  functions,  and  their 
first  derivatives  were  calculated  by  using  a  series  or  an 
asymptotic  expansion  depending  on  the  value  of  the  argument. 


-f  5 
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If  the  value  of  the  argument  was  in  the  range 

k 

0^rj^8  power  series  were  used  and  if  ^>8  the  asymptotic 

k 

expansions  were  used.  The  asymptotic  expansions  were 
needed  since  the  series  used  diverge  rapidly  for  arguments 
above  8.  The  values  of  the  Thomson  functions  given  by  the 
expansions  in  the  range  S^rj^lO  are  not  precise  and  vary 
up  to  approximately  10  per  cent  from  those  found  in  pub- 

*  k 

lished  tables  .  However,  this  error  was  not  noticeable 
in  the  portion  of  the  deflection  curves  calculated  using 
these  expansions. 

The  program,  which  was  run  on  the  University  of 
Alberta's  IBM  7040/1401  system  could  be  improved.  The 
series  and  asymptotic  expansions  used  were  computed  three 
different  times  in  the  program.  These  series  and  expansions 
should  have  been  in  a  subroutine  and  called  for  when  values 
were  desired.  A  recursion  formula  could  have  been  used 
in  place  of  actually  including  the  series  terms  in  the 
program.  This  would  save  computer  time  as  well  as  in¬ 
creasing  the  accuracy  of  the  results. 


*  Obtained  from  Handbook  of  Mathematical  Functions,  U.S. 
Department  of  Commerce,  National  Bureau  of  Standards, 
Applied  Mathematics  Series  55,  pp  384-385. 

**  Dwight,  H.B.,  Mathematical  Tables,  Dover,  1961,  pp  193- 
204. 
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An  explanation  of  the  symbols  used  for  the  input 
data  appear  below; 

BB  -  the  square  of  half  the  plate  width 

-6 

YM  -  Young's  Modulus  x  10 

NN  -  number  of  measuring  points  on  half  the  cross 
section  for  vh  ich  values  are  desired 
SMU  -  Poisson's  ratio 

TO  -  t 

o 

C  -  edge  thickness 
FLA  -  lambda  value 
OK  -  k 

The  results  appear  in  two  columns;  the  first  is 
the  non-dimensional  deflections  and  the  second  the  bend¬ 
ing  stresses  in  psi* 
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DIMENSION  Z(20) ,V( 30) ,B( 30) ,E( 30) ,D(20) 

COMMON  A { 6 , 6 ) » N  y  X ( 6 ) 

READ  (5,63)  B  B  ♦  YM 
63  FORMAT ( 2F7. 2 ) 

91  READ (5,92)  NN , SMU , TO , C , F L A , OK 

92  FORMAT  (  1 4 , 3F9 . 4 , 2F 7 . 3 ) 

WRITE (6,93)  FLA, OK 

93  FORMAT  ( 8H  LAMBDA=, F7  *3, 3  H  K=,F7.3) 

3  EP=2.*FLA*(C/(2.*TG)  ) *  *0 . 5 

ALP=2.*FLA*( 1 • +  C/ ( 2. *TQ )  )**0.5 

RAD=SQRT ( 3 . * ( 1.- SMU* SMU)  ) * 68/ ( FL A  * FL A *  TO ) 

8ET=C/ ( 2.*T0) 

P 1=3 . 14159 
R=l. 41421 
P=P I / 8 . 

EZ=EP/8 • 

AZ  =  AL  P/ 8  * 

IF(EP.GE.8. ) GO  TO  6 
DC  76  KK=2 ,28,2 
76  E ( KK ) =EZ* *KK 

ELOG=-ALGG( EP/2.  ) 

ABE=1.-64.*E(4)+113.77778*E(8)-32.36346*E( 12 ) +2 . 6419 1*E ( 16 ) 
BERE=-.08350*E(2Q )+.00123*E( 24)-.00001*E (28)+ A BE 
A8I=16.*E(2)-113.77778*E{6)+72.81778*£ ( 1 0 ) - 1 0 . 56 766*E (14) 
BEI£=.52186*£( 18 )- . 01 104*E ( 22 ) + .000 1 1 *E ( 26 ) +AB I 
ASK=EL0G*BERE+PI/4.*BEI£-.57722-59.0582*E (4) +171.36272*6 (8) 

BSK=— 60. 6097  7 *E (12)+5.6  55  39*£(16)-.  19 637* E ( 20 )  +  . 00 3 1 *E ( 24 ) 

SKEE=-. 00002 *E (28) +ASK+8SK 

DSK=ELOG*BEIE-PI /4.*BERE+6. 764 55* £ ( 2 ) - 1 42 . 9 1 828* E ( 6) 

ESK= 124. 23570* E( 10 ) -2 1.3006* E( 14 ) +  1 . 17  509*E ( 18 )-.02696*E (22 ) 
SKI£=.G0030*E(26  )+DSK  +  ESK 

BER1=EP* (-4. *6(2 ) + 14 . 22222* E ( 6 )-6 . 068 1 5* £ ( 10 ) +.66048*E ( 14) ) 

BRPE=6P* ( -.02609* E( 18)+.00C46*E(22))+B£R1 

BE  I 1=EP* (0.5-10.66  667*E(4 )  +  11.3  77  78*E { 8 )-2. 31 168*E ( 12 )  ) 

BIPE  =  £P* ( .14677*£( 16 )-.00379*E { 20  )  +  .00004*E ( 24  )  )+B£Il 
SKRP1=EL0G*BRPE-BERE/£P+PI/4.*BIPE+EP* ( - 3 . 69 1 1 4* E ( 2 ) ) 

SKRP2=EP*( 21.42034*E ( 6 )— 1 1. 36433* E ( 1 0 ) + 1 . 4 1 38 5*E ( 14) -. 06163 *E ( 1 8 ) ) 
SKRPE=EP* ( .001 16*E(22 ) - . 000 0 1  * E ( 2 6 )  )  +  SKRP  1  +  SKRP2 
SKI2=EL0G*BIPE-BEIE/EP~PI/4.*8RPE+.21 139*EP 

SKI3  =  EP*(-13.39859*E( 4 )  +  19 . 4 1 1 83* E ( 8 ) -4 . 6595 1 *£ (  12 )  +  . 33049 *E ( 16)  ) 
SKIPE  =  EP*(-.00927*E( 20 )  +  . 000 1 1 *E ( 24 ) ) + SK I  2+ SK I  3 
GO  TO  5 
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6  TPE=(2.*PI*EP)**(-0.5) 

PSE=SQRT (PI/ ( 2. *£P ) ) 

EXO=EXP ( EP/R ) 

EXH=EXP(-EP/R  ) 

BERE=TPE*EXO*CGS  (  EP/R-P) 

B£I£=TPE*EXO*SIN( EP/R-P  ) 

SKEE=PSE*EXM*CQS ( EP/R  +  P  ) 

SKIE=-PSE*EXM*SIN( EP/R  +  P  ) 

BRPE=TPE*EXG*CGS (EP/R  +  P  ) 

BIPE=TPE*EXC*SIN( EP/R+P) 

SKRPE=-PSE*EXH*CCS( EP/R-P) 

SKIPE=PSE*£XM*SIN(EP/R-P  ) 

5  I F( ALP.GE.8. ) GO  TO  8 
DC  74  J  J  =  2  »  28  »  2 
74  B  (JJ  )  = AZ*  * J  J 

PLGG=-ALGG( ALP/2.  ) 

ABA=1.-64.*B(4 )+ 1 13 . 77778*8 ( 8 )-32 . 3634 6* B ( 12) +2, 64191*6 ( 16) 
BERA=-. 08350*8(20  ) + .00 12 3*8 ( 24 )-. 0000 1*B ( 28  )+A6A 
ACI=16.*B< 2)-113.77778*B(6)+72.81778*B( 10 )- 10. 56766*8 ( 14) 
BEIA=.52186*8( 18 ) - . 0 1 1 04 * B ( 22 )  +  . 0 00 1 1  * B ( 26 ) + AC  I 
XSK=PLOG*BERA+PI /4.*BE I A- . 57722-59. 05820*B ( 4 ) 

YSK= 171.36272*B( 8 ) -60 . 60977 *B ( 1 2 ) + 5 . 65539*6 ( 16 )-. 19637*6 ( 20 ) 
SK£A=.0031*B( 24)-.00002*B(28 )+XSK+YSK 

ZSK=PLGG*BElA-PI/4.*B£RA+6. 7645 5* B(2)-142. 91828*6(6) 

ETK= 124. 2357*8 ( 1 0 ) -2 1 . 3006*8 ( 14 )  +  1 . 1 7 5 09*B ( 18 ) - . 02696*8 ( 22 ) 

SKI A=. 00030*8 ( 26  )+ZSK  +  ETK 

8EP1=ALP* l-4.*B( 2 )+ 14 . 22222*8 ( 6 )-6 . 063 15*B ( 1 0 ) + . 66048 *8 ( 14) ) 
BRPA=ALP* (-.02609*8 ( 1 8 ) + . 0C046*B ( 22 ) ) +BEP1 
BIP1= ALP *(0.5— 10.66667*8(4 )+ 11.37778*6 ( 8)-2.31168*B(12)) 
BIPA=ALP*(.14677*8(16 )-. 00379*8 ( 2 0 )+. 00004*6 ( 24 )  )+8IPl 
SPA=PLQG*BRPA-8ERA/ALP+P 1/4. *81 PA+ ALP* (-3.69114*8(2) ) 

TP=ALP*( 21.42034*8 (6)- 11. 36343 *B(  1 0 )  +  1 . 4 138 5* B {  1 4 ) - . 06 1 63* B {  18) ) 
SKRPA=ALP* ( .001 l6*B ( 22 )— .0C001*B( 26) )+SPA+TP 
EPA=PLGG*BIPA-BE I  A/ AL P-P I /4 . *BRPA  +  .21 1 39* ALP 

FP=ALP*(-13. 39859*6(4 ) +19.41183*6 ( 8 )-4 . 6595 1 *B ( 1 2 ) + . 33049*B ( 16) ) 
SKIPA=ALP*(-.00927*B( 20) +.C0011*B (24) J+EPA+FP 
GO  TO  9 

8  TPA=(2.*PI*ALP)**(— 0.5) 

PS A= SORT ( PI/ ( 2.*ALP  )  ) 

E XA=E XP ( ALP/R  ) 

EXB=EXP{— ALP/R  ) 

BERA=TPA*EXA*CGS ( ALP/R-P  ) 

BE  I A=TPA*EXA*SIN( ALP/R-P ) 

SKEA=PSA*EX8*CGS ( ALP/R +P ) 

SKI A=-PSA*EXB*SIN( ALP/R+P) 

BRPA=TPA*EXA*COS (ALP/R+P  ) 

B IPA=TP A *EXA* SIN (ALP/R+P ) 

SKRPA=-PSA*EXB*COS( ALP/R-P) 

SK I PA=PSA*£XB*SIN( ALP/R-P) 
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9  A  ( 1  , 1 )=-ALP*BEIA-2.*BRPA 
A { 1,2 )=ALP*B ERA-2. *BIPA 
A ( 1 ,3 )=-ALP*SKIA-2.*SKRPA 
A ( 1 , A  )  =  ALP*SKEA-2.*SKIPA 

A ( 1 ,5)=+2.*T0*ALP**3/ (4.*FLA*FLA) *( l.+OK) 

EP2=EP*EP 
EP3=EP**3 
A  LP2  =  ALP  * ALP 
ALP3=ALP**3 

A (2, 1  )=4.*EP*BEIE-EP2*BIPE+8.*8RPE 
A(2,2)*-4.*EP*BERE+EP2*BRPE+8.*B  I  PE 
A(2,3)=4.*EP*SKIE— EP2*SKIP£+8.*SKRPE 
A( 2,4)=-4. *EP*SKEE+EP2*SKRPE+8.*SKIPE 
A(2»5)=+SMU*BB*BET/(RAD*FLA*FLA)*EP3 
A( 3, 1 )=-EP3*BERE-24.*EP*BE IE-48. *88PE+8.*£P2*BIPE 
A ( 3 ,2)*-EP3*BEIE+24. *EP*BERE-48. *8IPE-8.*EP2*BRP£ 

A ( 3 , 3 )=-EP3*SKEE-24. *£P*SK IE-48. *SKRPE+8.*£ P2* SKI  PE 
A(3,4)*-EP3*SKIE+24.*EP*SKEE-48.*SKIPE-8.*EP2*SKRPE 
A  (  3 , 5  )  =  0  . 

A (4, 1)=ALP2*B ERA-2. *ALP*B  I  PA-EP2*8ERE+2.*EP*8I  PE 
A(4,2)=ALP2*BEIA+2.*ALP*8RPA— £P2*B£IE— 2.*£P*BRPE 
A (4, 3)=ALP2*SKEA-2.*ALP*SKIPA-£P2*$KEE+2.*EP*SKIP£ 

A(4,4)=ALP2*SKI A+2. *ALP*SKRPA-EP2*SKIE-2.*EP*SKRPE 
A(4,5)= (-1.5) *SMU*TG*( ALP**4-EP**4)/ ( FLA*FLA*SQRT ( 3.-3. *SMU*SMU ) ) 
N=4 

CALL  SOLVE 
NW=NN+1 
DO  99  J= 1 , NW 
BJ=J— 1 

BN*NN 

ETA*2.*FLA*SGRT(  1 . -B J / 8N+C / ( 2  .  *TO  )  ) 

V  A=ET  A/8 . 

I F { ETA.GE.8 . ) GO  TO  12 
DO  78  K=2,28,2 
78  V { K )  =  VA*  *K 

XLOG=— ALOG (ETA/2 . ) 

BPN1*ETA*(-4.*V(2)+14.22222*V(6)-6.06815*V( 10 )+. 66048* V ( 14)  ) 
BRPN=ETA*(~. 02609* V( 1 8 ) + . C 0 046*V ( 22 ) )+BPNl 
BPN2=ETA*(0. 5— 10.66667*V(4)+11.37778*V(8)— 2.31168*V(12) ) 

B  I PN=  ET A* ( . 14677 *V (  16  ) - . 00 3 79 * V ( 2 0 ) + . C0004* V ( 24 )  )+BPN2 
BRA=1. -64. *V(4)+113.77778*V(8)-32. 36346* V(12)+2.64191*V(  16) 

BERN*-. 08350* V( 20 ) +  .00  123* V ( 24) -.00001* V (28 )+BRA 
B IN* 16 . * V ( 2  )-l 13. 7777 8* V( 6 )+72.81 778*V ( 10 )- 10. 56  766*V (14) 
BEIN*.52185*V(18)-.01104*V( 22)+.00011*V(26)+BIN 
S PN= X LOG * B RPN— BERN/ ET  A  +  PI/4.*BIPN+ETA* ( - 3 . 69 1 1 4* V ( 2 ) ) 

TNT*ETA* (21.4  2034*V( 6  )-ll. 36343 *V(10)+1. 41385* V(14)-. 06 163* V(18)) 
SKRPN  =  ET A* (  .0Q116*V ( 22  )-.0CC01*V( 26) ) +  SPN  +  TNT 
EPN=XLOG*BI PN—  BE  IN/ETA-PI/4. *BRPN+ETA*. 21139 

FN=ETA* (~13.39859*V (4 )  +  l9.41183*V { 8 ) -4 . 6595 1  * V ( 1 2 )  +  . 3 3049* V (16)) 
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SKIPN=ETA*(-.G0927*V{ 20 ) + . C CO 1 1 *V ( 24 ) ) +EPN+FN 
XSK=XL0G*BERN+PI/4.*BE  IN- . 577  2  2— 59 .05  8  2* V ( A )  +  1 71  ,36272*V(  8) 
YS=-60.60977*V (  12)+5.65539»V( 1 6 )- .  1 9637* V ( 2 0 ) + . 003 1 *  V ( 24) 
SKEN=— ,0GQ02*V (28 )+XSK+YS 

DEK=XL0G*BEIN-PI/4.*BERN+6.76455*V(2)-142.91828*V(6) 

EEK= 124, 23570# V( 1 0  ) -2 1 . 3006* V { 1 4 )  +  1 . 1 7 509* V { 1 8 ) - . 02696 * V ( 22 ) 
SKIN=.00030*V (26 )+DEK+EEK 
GC  TO  13 

12  TPN=(2.*PI*ETA)**(-0.5) 

PSN=SQRT ( PI/ ( 2.*ETA ) ) 

EXC=EXP ( ETA/R ) 

EX D=EXP(~ ETA/R ) 

BRPN=TPN*EXC*CGS ( ETA/R+P ) 

BIPN=TPN*EXC*SIN (ETA/R+P ) 

SKRPN=-PSN*EXD*CQS (ETA/R-P } 

SKIPN=PSN*£XC*SII\i( ETA/R-P ) 

BERN=TPN*EXC*CGS (ETA/R-P ) 

BE  I N-TPN*EXC*S I N ( ETA/R-P ) 

SKEH=PSN*EXD*CGS( ETA/R+P ) 

SKIN=-PSN*EXD*S INI ETA/R+P ) 

13  F=C/ ( 2 • *  TO ) 

G AMK=  TO# ( 1  •  — BJ/BN  +  F ) * ( 1 , +0K } -OK *C /2. 

CGNT=TQ* I  2  *  * (  l.+0K)/3.  +  l.*F*( 2 . +0K ) +2 . *F*F ) / ( 2 . * F+ 1 . ) 

WC=-1.* ( GAMK-CCNT ) 

VLT=6.*SMU*TG/(FLA*FLA*SQRT(3.-3, *SMU*SMU) ) 

Z(J)  =  1./ETA*(X(  1  )*BRPN  +  X(  2  )*8IPN+X(  3)  *  SKRPfsi  +  X  (  4  )  *SKIPN)-W0-VLT 
D(J)=(Z(J)+WG)/(T0+C/2. ) 

DIMENSION  W ( 5 ) , S X ( 2 0 ) 

T  =  2  . *T0* ( 1 BJ/BN ) +C 
DG  67  KE  = 1 1  5 
67  W ( KE ) =ETA**KE 

F  TD=X ( 1 ) * ( 4 , * W ( 1  )*BEIN-W(2)*8IPN+8.*BRPN) 

STD=X (2) * (-4,*W( 1 }*BERN+W( 2 ) *BRPN+8,*B IPN ) 

T  T  D=  X ( 3 ) * ( 4 • *W ( 1) *SKIN— W(2)*SKIPN+8,*SKRPN) 

ATD=X (4 ) * (-4,*W( 1 ) *SKEN  +  W ( 2 ) *SKRPN+8,  *SKI PN ) 

CGN=4.*FLA**4/ ( BB*W ( 5 ) ) 

SDZ=CCN* (FTD+STD+TTD+ATD) 

FC=+YM* ( 1 ,E06 ) *T / ( 2 • * ( l.-SMU*SMU> ) 

SX( J )=FC* { SDZ+SMU/R AD } 

WRITE  (6,45)  D ( J ) »  SX ( J ) 

45  FGRMAT( lXfF10,5, 10XtF10.2) 

99  CGNTINUE 
GO  TO  91 
END 
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SUBROUTINE  SOLVE 
COMMON  A (6,6)  ,  N , X ( 6 ) 

LOGICAL  SING 
S  ING=. FALSE. 

M=N+ 1 
N  M 1  =  N  - 1 
DC  I  1=1, N 

I  A ( I »  M ) =- A { I , M  ) 

DC  20  1=1, NM1 
AMAX  =  A ( I , I  ) 

K=  I 

DC  19  J=  I  ,  N 

IF  (ABS  (  A{  J  ,  I  )  )  .LE.ABS  UMAX  ))  GO  TO  19 
AMAX  =  A( J, I ) 

K  =  J 

19  CONTINUE 

IF  {  AMAX. NE. 0.0)  GO  TO  30 
SING  =  .FALSE. 

RETURN 

30  IF  (  ASS ( AMAX ) .GE .0. 16-20 )  GO  TO  31 
WRITE  (6,14) 

14  FORMAT  ( 1HJ , 14X , 45H  MATRIX  IS  EITHER  SINGULAR  OR  ILL  CONDITIONED) 

31  DO  18  J=i,M 
TS  =  A  (  I  ,  J  ) 

A  l  I  ,  J  )  =  A  {  K  ,  J  ) 

18  A ( K , J )  =  T  S 

IP1  =  1+1 
DO  11  J=  I  P 1 , N 
FMULTP  =  A( J, I ) / A { 1,1) 

DO  11  L=  I P 1  *  M 

II  A ( J , L  )  =  A ( J , L  )  -  FMULTP  *  A  <  I , L ) 

20  CONTINUE 

X (N )  =  A(N,M)/A(N,N) 

DO  22  1=1, NM1 
K  =  N-I 
KP1  =  K+l 
FNUM  =  0.0 
DO  23  J  =  K P 1 , N 

23  FNUM  =  FNUM  +  X(J)*A(K,J) 

22  X(K)=(A(K,M)-FNUM)/A(K,K) 

RETURN 

END 
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APPENDIX  II 

FORTRAN  IV  PROGRAM  FOR  CALCULATION  OF  EXPERIMENTAL  RESULTS 

The  IBM  7040/1401  was  again  used  to  numerically 
integrate  the  curvatures  obtained  from  tests  on  the  two 
bi-trapezoidal  plates.  The  procedure  is  outlined  in 
Chapter  II. 

An  explanation  of  the  symbols  used  for  the  input 
data  appears  below: 

NT  -  number  of  measuring  stations  on  half  the 
cross  section. 

T ( 1Z) -  thickness  of  the  plate  at  measuring  station 
IZ 

E  -  edge  thickness  of  the  plate 
NI  -  NT  -  1 

TT(JK)  -  thickness  of  the  plate  midway  between 
measuring  station  JK  and  JK  +  1 . 

YM  -  Young ' s  Modulus  x  10  ^ 

SMU  -  Poisson's  ratio 

HTT  -  half  maximum  plate  thickness 

BB  -  the  square  of  half  the  plate  width 

OK  -  k 

DG  -  distance  between  measuring  stations 

RF  -  range  factor  of  the  digital  data  processor 

(1.0  when  the  strain  is  in  micro  inches  per  inch) 
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I(K)  -  initial  strains  read  with  the  plate  on  a 
flat  surface 

J(LK)-  strains  recorded  under  load 

CC  -  E 

FLA  -  lambda  value. 

The  strains  on  one  half  of  the  cross  section  were 
recorded;  first  the  upper  side  starting  from  the  plate 
center  then  the  lower  side  again  starting  from  the  center. 
The  strains  on  the  other  half  were  then  recorded  in  a 
like  manner.  Since  there  were  26  gauges  on  half  the 
cross  section ,  gauges  1  and  14  were  on  the  top  and  bottom 
sides,  respectively,  at  the  cross  section  center.  These 
number  1  and  14  gauges  were  re-read  in  recording  the  sec¬ 
ond  half  of  the  cross  section. 
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DIMENSION  G(30),T{15),C(20),S(20),W(20),WT(20),WC(20),F(20),TT(20 
DIMENSION  CA(20),EXB{20),SXB(20) 

DIMENSION  1(30), J(30) 

READ (5,4)  NT  ,  T  0 , (T(  IZ),IZ=1,NT) 

4  FORMAT  { 15, F8.5/1 1QF7.4) ) 

RE AD ( 5 »  7 )  E,NI, (TT( JK), JK=1,NI) 

7  FORMAT  ( F8.5, 15/ {  10F7.4)  ) 

READ (5, 5)  YH,SMU,F)TT,BB,OK,DG,RF 
5  FORMAT  ( 5F7.3,F7.5,F7.3) 

2  RE AD ( 5  »  55 )  ( I ( K  )  , K= 1,26) 

55  FORMAT  (  ( 1315)  ) 

READ (5,56)  ( J ( LK ) , LK=  1 , 26 ) 

56  FORMAT  <  (1315)  ) 

READ (5*47)  CC,FLA 

47  FORMAT {2F7. 3) 

DO  57  KM=  1,26 
RI= 1 1  KM ) 

RG=  J ( KM ) 

57  G ( KM ) = ( RG— R I ) *RF 

RAD=SQRT ( 3, * ( l.-SMU*SMU) )  *BB/ ( FLA *FLA* TO ) 

DC  8  11=1,13 

EX6( I  I )={G( I  I )-G( I  1  +  13  )  )/2. 

SXB ( I  I )  =  YM/ { 1 SMU*SMU ) * ( EXB ( 1 1 )  +  SMU* 1 . 0E06*T C II  )/{2.*RAD)  ) 

8  CIII )  =  (G( II )-G( I  1  +  13) )/T(II ) 

S { 1 )=0. 

DO  9  L=2 , 13 

9  S (L)  =  S(L-1 )+0.5*DG*(C ( L ) +C ( L-l)  ) 

W ( 1 )=G. 

DC  11  M=  2 , 1 3 

11  W(M)=WIM-1)+.5*DG*(S(M)+S(M-1) ) 

SCA=0.0 

A=SORT ( BB ) * ( TO+E ) 

NB=N I - 1 
DC  19  N J=2 , NB 

CA(NJ)=TT(NJ)*(W(NJ)+W(NJ+1)  ) 

19  SCA=SCA+CA(NJ ) 
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WBA=0.5*0G*<TT ( 1  )* (W(  1)+W(2) )+SCA+TT(NI )*<W(NI)+W<NI+1))) 
W8=W8A/A 
DC  13  IA=1,13 
13  WT ( I  A  )  =  W  ( I A )  —  WB 
B=CC/ ( 2 . *T0  ) 

DO  15  NN=0  t 12 
B  J  =  NN 
BN=12.0 

GAMK= ( l.-BJ/BN  +  B  )*TG* ( l.+OK )-OK*CC/2. 

CGNT=TO* (2. * ( l.+OK  ) /3.+l.*E* ( 2.+OK ) +2 . *B*8 ) /{ 2.*B+1. ) 

15  WC(NN+1)=(GAMK-CQNT)*{-1.0) 

DO  17  N A= 1 ,  1 3 

F(NA)=(  (1.0E-Q6)*WT(NA)+W0<NA)  )/HTT 
17  CONTINUE 

WRITE  (6,99)  FLA, (F{ JJ  )  , JJ=1, 13) 

99  FORMAT ( 8H  L AMBDA= , F7 . 3/ ( IX , 5F IQ. 6 ) ) 

WRITE  (6,98)  (SXB(LZ  )  , LZ=  1 »  13  ) 

98  FORMAT ( 17H  BENDING  STRESSES/ (  IX , 5F 12. 2  )  ) 

GO  TO  2 
END 
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APPENDIX  III 

DETAILS  OF  CARD  PUNCH  USE 

An  IBM  024  key  punch  (card  punch)  was  used  to  re¬ 
cord  the  strains  which  were  measured  by  the  digital  data 
processor.  Figure  A3-1  shows  the  modifications  and  elect¬ 
rical  connections  to  the  key  punch  necessary  for  operation 
with  the  processor.  The  connections  shown  were  wired  to 
a  "Cannon"  connector  which  was  plugged  into  the  key  punch 
outlet  in  the  processor.  The  numbers  shown  in  figure 
A3-1  refer  to  those  on  the  terminal  board  of  the  punch 
and  are  shown  again  on  IBM  Wiring  Diagram  No.  228001P. 

In  order  for  the  cards  punched  to  be  compatible 
with  the  program  of  Appendix  II  a  drum  card  for  the  IBM  024 
was  used.  This  card  is  shown  in  figure  A3 -2.  With  this 
card  thirteen  channels  of  information  were  recorded  per 
data  output  card. 


73 


CANNON 


FIGURE  A3-1  KEY  PUNCH  TERMINAL  PANEL  MODIFICATION 
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FIGURE 
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APPENDIX  IV 

MEASUREMENT  OF  MOMENT  VERSUS  CURVATURE  RELATION 

The  moments  required  to  produce  the  curvatures 
were  measured  and  recorded.  These  are  shown  in  figures 
A4-1  and  A4-2.  Due  to  eccentricities  in  the  loading 
apparatus,  the  moments  recorded  are  only  approximate.  The 
moment  needed  to  counterbalance  the  weight  of  the  plate 
to  bring  the  curvature  initially  to  zero  has  been  sub¬ 
tracted  from  all  values. 
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FIGURE  A4-1.  MOMENT  VERSUS  A  FOR  PLATE 
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FIGURE  A4- 2 .  MOMENT  VERSUS  A  FOR  PLATE  II 


